Reconciling the Light Scalar Mesons 
with Breit-Wigner Resonances as well as the Quark Model 



Eef van Beveren 

Centra de Fisica Tedrica 
Departamento de Fisica, Universidade de Coimbra 
P-3000 Coimbra, Portugal 
eef @teor .fis.uc.pt 

George Rupp 

Centro de Fisica das Interacgoes Fundamentais 
Instituto Superior Tecnico, Edificio Ciencia 
P-IO49-OOI Lisboa Codex, Portugal 
george@ajax.ist.utl.pt 

PACS number(s): ll.80.Et, 12.40.Yx, 13.75.Lb, 14.40.-n 

February 1, 2008 



Abstract 

Resonances appearing in hadronic scattering processes are described by a two-phase 
model. In the one phase, scattering products are observed, whereas the other phase de- 
scribes confinement. A so-called "Resonance-Spectrum Expansion" is derived, containing 
expressions that resemble Breit-Wigner formulae. This method also provides a straightfor- 
ward explanation for the origin of the light scalar mesons without requiring extra degrees 
of freedom. 



1 Introduction 

For more than three decades now, the light scalar mesons have been puzzling both experimen- 
talists [1-4] and theorists [5-12]. On the experimental side, one is faced with a highly disparate 
set of a few well-established, relatively narrow resonances, as well as some very broad, seemingly 
non-Breit-Wigner-like structures in S'-wave meson-meson scattering. On the other hand, theory 
appears to have tremendous difficulties in reproducing these states on the basis of some micro- 
scopic quark substructure. For instance, the naive quark model would describe these states as 
confined P-wave qq systems, thus putting the masses of the lightest nonet at about 1.3 GeV 
upwards, and not below 1 GeV as found experimentally. 
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The first consistent attempt to describe the lowest scalars as a usual meson nonet was un- 
dertaken by R. L. Jaffe in 1977, in the framework of the MIT bag model [13]. In this approach, 
the lightest exotic q^q^ scalar states indeed form a nonet, and are therefore also called crypto- 
exotics. Moreover, due to a very large, attractive color-magnetic interaction term, the central 
mass values of the corresponding bag states are shifted downwards several hundreds of MeV, 
thus being in rough agreement with the real parts of the light scalar mesons listed in the PDG 
tables [14], including the /o(600) or a meson, in those days called the e. However, it is not at 
all clear how to couple these stable multiquark states to the physical thresholds, and especially 
what influence this would have not only on the widths but also on the real part of the spectrum. 
Then, in 1982, M. D. Scadron presented the first theoretical description of the light scalars as 
a nonet of qq states [15]. In this work, the mechanism for producing light scalar mesons is the 
spontaneous breaking of chiral symmetry, at the same time responsible for the vanishing of the 
pion mass in the chiral limit. However, with the advent of chiral perturbation theory (ChPT), 
the (T and k mesons fell out of favor, and for many years it was considered unnecessary, perhaps 
even undesirable, to have a complete light scalar nonet. Only with the undeniable mounting of 
experimental evidence, first for the /o(600) [16] and recently also for the Kq{^Q{)) (k) [17], many 
model builders started to rehabilitate the light scalar nonet, even some ChPT practitioners [18]. 
This development also led to a revival of interest in our 1986 model prediction of the complete 
light nonet, in a unitarized Schrodinger formalism [19]. 

In modern approaches towards fundamental interactions and scattering of elementary parti- 
cles, the Schrodinger equation has become obsolete. However, as we may learn from the long- 
standing difficulties involving the light scalar mesons, an analysis based on the Schrodinger 
equation, like the one employed in Refs. [19, 20] can be very clarifying. Unfortunately, it is 
then necessary to enter into the details of solving the corresponding set of coupled second-order 
differential equations, a technique which has become obsolete as well. 

In Table 1, we give a simple classification for the positive-parity mesons, based on the model 
we describe in Sec. 2. Each resonance is supposed to originate from a pure confined qq state, 
with quantum numbers s = 1 for the total intrinsic qq spin, and £ = 1 or 3 for the relative orbital 
angular momentum of the qq pair. In Sec. 2 we explain what we exactly mean by "originating 
from", and how this can be consistently described within a model for meson-meson scattering. 
By just varying one model parameter, i.e., the coupling A, letting it decrease from its "physical" 
value Xph towards zero, the model's spectrum turns from the experimentally observed resonances 
into a genuine qq confinement spectrum. 

For the radial confinement spectrum we have chosen equal spacings in Table 1, with level 
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Table 1: The classification of the = (0, 1,2)+ resonances as a function of radial excitation n. The 
ground states have n = 0. 

splittings of 350 MeV and a ground state {n — 0) at 1.4 GeV. The i — 3 radial excitations start 
out at 1.75 GeV. Such a spectrum corresponds to the one of a harmonically oscillating quark- 
antiquark pair. Harmonic confinement is not essential to our model, but the cross sections, 
phase shifts and electromagnetic transition rates resulting from the full model are in reasonable 
agreement with experiment [19,23,24]. 

Non-strange (nn) and strange (ss) configurations double the number of isoscalar states into 
SU (3)-fiavor singlets and octets. But one should be be aware that all states get mixed through 
the meson loops in our model (see Sec. 3). Hence, like in Nature we will not obtain pure angular, 
radial, or fiavor excitations in our model calculations. 

In Table 2 we have schematically indicated how many states one must expect at each mass 
level, except for the ground states where i — 3 excitations are absent. 

Meson loops infiuence, moreover, the precise resonance shapes. Some come out broad, others 
narrower. Also, the central resonance positions may shift substantially (100-300 MeV) with 
respect to the underlying qq confinement spectrum. 

Prom Tables 1 and 2 we may conclude that the observed positive-parity mesonic resonances 
can easily be accommodated in a quark model, contrary to what has been claimed in recent 
literature [25,26]. Notice in particular the claim in Ref. [25] that "far too many resonances 
were established, to be accommodated in the ground-state scalar nonet", a few lines further on 
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Table 2: The expected number of = (0, 1, 2)+ meson resonances at each level of radial excitation. 

followed by the remark "... a missing state". Inspection of Table 1 reveals that still many states 
are missing, especially in the 1 = 1 and 1 = 1/2 sectors. However, thanks to glueball searches 
we nowadays have a much better knowledge of the 7 = sector than twenty years ago. A 
classification of the mesonic resonances in this sector can be satisfactorily achieved assuming 
quark degrees of freedom. That does not necessarily imply the absence of other configurations. 
From lattice calculations it is becoming clear that it may even be very hard to disentangle the 
various configurations existing in /q systems [27]. Hence, the /qS could, in principle, be mixtures 
of qq, glueballs [28], {qq)^, meson-meson states, hybrids, etc., moreover in all possible color 
configurations. Nevertheless, for the classification of positive-parity mesonic resonances we only 
need quark degrees of freedom. 

However, our table does not contain the light scalar mesons, that is, the /o(600) {a meson), 
the /o(980) and Oo(980), and the recently confirmed [17] irQ(800) {k, meson). But we have 
shown that such a complete nonet below 1 GeV is inevitable in the quark model, due to the 
^Po mechanism [19,20], besides the usual nonet in the mass region 1.3-1.5 GeV, consisting of 
the /o(1370), ao(1450), K*{U30), and /o(1500). Notwithstanding the dissent about the correct 
interpretation of the individual states, these two scalar nonets arc the only ones that seem to be 
complete from the experimental point of view. Several members of the other nonets still have 
to be found. Table 1 shows in which energy region we expect the missing qqS and KqS to be 
observed in experiment. 

Let us now dwell somewhat more upon the understanding of the scalar mesons, and compare 
their situation with positronium. Imagine a teacher asking his/her students, after having added 
a fictitious ground state, to invent a theory that explains the positronium spectrum. Probably 
only the cleverest students will discover the mahce of their teacher. The others will invent 
whatever model it takes to get an explanation for the false positronium spectrum. In the scalar- 
meson spectrum, at least two such extra states seem to exist, namely the firmly established, 
relatively narrow /o(980) and Oo(980) resonances, which are therefore also the most controversial 
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states for theorists. As an illustration of the general confusion here, let us just mention e.g. 
the "XX-molecule" approach of Ref. [29], producing both the /o(980) and ao(980) as dynamical 
meson-meson resonances due to strong t-channel attraction in the KK system, the relativistic 
quark model of Ref. [30], in which only the /o(980) is described as a qq state, owing its low mass 
to a strong instanton-induced interaction, whereas the ao(980) is supposed to be a KK molecule, 
and finally the confining NJL model of Refs. [31,32], which obtains the /o(980) as well as the 
ao(980) as light qq scalars, thanks to an attractive 't Hooft interaction (sec, however, Ref. [33]). 
Without these two resonances, the scalar meson spectrum would seem to start off with a ground 
state at about 1.3 GeV and could then be simply explained by the naive quark model, provided 
one assumed the very broad (7(600) and «;(800) resonances to be of whatever dynamical, but not 
qq origin. On the other hand, if one insists on taking the light scalars as the ground states of 
the spectrum, then, at first sight, it seems very hard to achieve a unified description of all scalar 
mesons, not to speak of including pseudoscalar, vector, and other mesons, too. 

The main purpose of the present paper is to demonstrate that it is really possible to uniformly 
describe the whole nonet of light scalar mesons, including the broad a and n structures, in a ^Pq- 
modified Breit-Wigner-like framework, on the basis of qq and meson-meson degrees of freedom 
only. However, these states will turn out to be not the naive ground states of the scalar-meson 
spectrum [19,20,34]. The formalism we shall use is just the old-fashioned Schrodinger equation, 
often considered "unworthy" in hadronic physics nowadays, but which nonetheless is perfectly 
suited to describe phenomena like (virtual) bound states, resonances, threshold behavior, Rie- 
mann sheets, and so forth. Of course, we do not propose this technique so as to substitute the 
modern nonperturbative methods for handling QCD. It is just intended to pinpoint the structure 
of the scalar-meson spectrum, and to show how powerful a simple, intuitive approach can be, as 
long as one includes the relevant degrees of freedom. For that purpose, one unavoidably has to 
go through some widely forgotten calculus, in order to obtain an analytic expression for the S 
matrix from a set of coupled-channel equations. Therefore, we shall choose the simplest possible 
interactions that still contain the essence of the physics, without eclipsing the latter by a heap 
of opaque equations. 

The organization of this paper is as follows. In Sec. 2 we formally present the coupled- 
channel equations linking the confined qq sector to the free two-meson sector. Section 3, in 
combination with Appendices A, B, and C, is devoted to a detailed analytical derivation of the 
S matrix and the scattering phase shift for a two-channel model, with arbitrary confinement and 
a delta-shell interaction to mimic the ^Pq transitions between the qq and meson-meson channels. 
The here derived formula, which we call "Resonance-Spectrum Expansion", is central to the 
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ensuing analysis in the remainder of the paper. In Sec. 4 we discuss how to extract results for 
bound states, resonances, and scattering observables from the S matrix for the full, multichannel 
model. In Sec. 5 we describe the general connection between the discrete confinement spectrum 
and the ^'-matrix poles corresponding to physical bound states and resonances, depending on 
the threshold energies. Section 6 explains the difference in threshold behavior of 5'-wave and 
P-wave poles. In Sec. 7 we present a detailed discussion of the behavior of 5'-matrix poles, as a 
function of the relative momentum k or the energy E, for the light scalar mesons, in particular 
the Kq{800) and the ao(980). Our conclusions on the nature of the scalar mesons are drawn in 
Sec. 8. 

2 Coupled channels 

Let us consider a mesonic system which may appear in two different phases, and a mechanism 
allowing transitions from one phase to the other. In one phase the system consists of two 
noninteracting mesons, in the other phase of two permanently bound particles representing a 
quark-antiquark system. We will refer to the former phase as free, to the latter as confined. 
The communication between the two phases we describe through a short-range potential Vt. In 
the interaction region, which has the dimension of about 1 fm, we allow both phases to coexist. 
Hence, the wave function of such a system consists of two components there, ipf and ipc- Outside 
the interaction region the confined component ipc must, of course, vanish rapidly. This can be 
achieved by a potential Vc, rising to infinity with distance, like the linear and the harmonic- 
oscillator potentials. 

Let the Hamiltonians He and Hf describe the dynamics of the phases of two permanently 
confined particles and of two free particles, respectively. Then the following set of coupled 
Schrodinger equations describes the dynamics of the full system: 

(E-H,) V^e(r) = Vt^fif) and {E - Hf) t/jfif) = [Vtf ^c{r) . (1) 

Such a philosophy already underlied an elaborate coupled-channel quark model [19,23,35,36], 
designed to simultaneously describe mesonic bound-state spectra, resonances, and meson-meson 
scattering. However, in spite of the model's success to reproduce a host of experimental data 
with a very limited number of parameters, it is not very suited for the point we wish to make in 
the present paper, owing to the specific model choice of the confining qq potential, and especially 
the rather complicated matrix expressions needed to obtain jS-matrix-related observables. Thus, 
we shall use here an arbitrary confinement potential, and a very simple transition potential Vf. 
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3 Quark exchange and meson loops 



When, as depicted in Fig .1, two interacting mesons exchange a quark, the resulting system will 
consist of a valence quark-antiquark pair. Whether this qq pair is going to form a resonance or 
not will depend on the quantum numbers of the system and the total available energy. 




Figure 1: The mesons (M) exchange a quark or, equivalently, a quark-antiquark pair is annihilated, 
followed by a second quark exchange, equivalent to a new quark-antiquark pair being created. 



Near a resonance, such a process may be described by scattering phase shifts of the form 



cotg{5i{s)) 



E 



R 



Tr/2 



(2) 



where Eji and F^^ represent the central invariant meson-meson mass and the resonance width, 
respectively. 

However, formula (2) is a good approximation for the scattering cross section only when the 
resonance shape is not very much distorted and the width of the resonance is small. Moreover, 
the intermediate state in such a process is essentially a constituent qq configuration that belongs 
to a confinement spectrum (also referred to as bare or intrinsic states), and hence may resonate in 
one of the eigenstates. This implies that the colliding mesons scatter off the whole qq confinement 
spectrum of radial, and possibly also angular excitations, not just off one single state [37]. Con- 
sequently, a full expression for the phase shifts of Eq. (2) should contain all possible eigenstates 
of such a spectrum, as long as quantum numbers are respected. Let us denote the eigenvalues 
of the relevant part of the spectrum by En {n = 0, 1, 2, . . .), and the corresponding eigenstates 
by Fn. Then, following the procedure outlined in Appendices (A), (B), (C), and in Ref. [20], we 
may write for the partial-wave phase shifts the more general expression 



cotg((5(s)) 



n=0 V * 



n -1 r 



Er, 



n=0 V-S 



T I 



En 



(3) 



which we call "Resonance-Spectrum Expansion". In R{s) and I{s) we have absorbed the kine- 
matical factors and the details of two-meson scattering, and moreover the three-meson vertices. 
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For an approximate description of a specific resonance, and in the rather hypothetical case 
that the three-meson vertices have small coupling constants, one may single out, from the sum 
over all confinement states, one particular state (say number N), the eigenvalue of which is 
nearest to the invariant meson-meson mass close to the resonance. Then, for total invariant 
meson-meson masses y/s in the vicinity of Ej^, one finds the approximation 



'En + R{s) 




2 




I{s) 


Tn 


2 



cotg{S{s)) « ^ ■ (4) 



Formula (4) is indeed of the Breit-Wigner form (2), with the central resonance position and width 
given by 

Er ^ En + R{s) \J^Nf and Tr ^ 2/(s) {J^^f ■ (5) 

In experiment, one observes the influence of the nearest bound state of the confinement spectrum, 
as in classical resonating systems. Nevertheless, Eq. (4) is only a good approximation if the three- 
meson couplings are small. Since the coupling of the meson-meson system to quark exchange is 
strong, the influence of the higher- and lower-lying excitations is not negligible. 

In the other hypothetical limit, namely of very large couplings, we obtain for the phase shift 
the expression 

cotg(5(.)) ^ ^ , (6) 

which describes scattering off an infinitely hard cavity. 

The physical values of the couplings come out somewhere in between the two limiting cases. 
Most resonances and bound states can be classified as stemming from a specific confinement 
state [38,39]. However, some structures in the scattering cross section stem from the cavity 
which is formed by quark exchange or pair creation [20]. The most notable of such states are the 
low-lying resonances observed in 5- wave pscudoscalar-pscudoscalar scattering [17,40]. 

From the discussion above one may conclude that, to lowest order, the mass of a meson follows 
from the quark- antiquark confinement spectrum. It is, however, well-known that higher-order 
contributions to the meson propagator, in particular those from meson loops as depicted in Fig. 2, 
cannot be neglected. 

Virtual meson loops give a correction to the meson mass, whereas decay channels contribute 
to the strong width of the meson, too. One obtains for the propagator of a meson the form 

n(5) = —— ^ , (7) 

* (-^confinement ^ ^-^meson loops j 




meson loop 



Figure 2: The lowest-order self-energy graph for a meson propagator. 

where AM develops complex values whenever the threshold of a decay channel gets surpassed. 

For the full mass of a meson, all possible meson-meson loops have to be considered. A 
model for meson-meson scattering must therefore include all possible inelastic channels as well. 
Although in principle this could be done, in practice it is not manageable, unless a scheme 
exists dealing with all vertices and their relative intensities (see e.g. Ref. [41] for scalar and 
pseudo-scalar meson couplings). In Refs. [42,43] relative couplings have been determined in the 
harmonic-oscillator approximation assuming '^Pq quark exchange. However, further kinematical 
factors must be worked out and included. 

4 The spectrum 

The full model consists of an expression for the K matrix similar to Eq. (3), but extended to 
many meson-meson scattering channels, several constituent quark- ant iquark channels, and more 
complicated transition potentials [19,23], which at the same time and with the same set of four 
parameters reproduces bound states, partial-wave scattering quantities, and the electromagnetic 
transitions of cc and bb systems [24]. 

The K matrix can be analytically continued below the various thresholds, even the lowest 
one, with no need to redefine any of the functions involved, in order to study the singularities of 
the corresponding scattering matrix. Below the lowest threshold, these poles show up on the real 
^/s axis, and can be interpreted as the bound states of the coupled system, to be identified with 
stable mesons. For the light flavors one flnds this way a nonet of light pseudoscalars, i.e., the 
pion, kaon, eta, and eta'. For the heavy flavors, the lowest-lying model poles can be identifled 
with the £1(1870), £,(1970), r]c{lS), J/^{IS), ^^(3686), 5(5280), 5,(5380), T(15), T(25), and 
T(35). 

Above the lowest threshold, the model's partial-wave cross sections and phase shifts for all 
included meson-meson channels can be calculated and compared to experiment, as well as the 
inelastic transitions. Singularities of the scattering matrix come out with negative imaginary part 
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in the ^/s plane. To say it more precisely: out of the many singularities in a rather complex set 
of Riemann sheets, some come close enough to the physical real axis to be noticed in the partial- 
wave phase shifts and cross sections. In fact, each meson-meson channel doubles the number of 
Riemann sheets, hence the number of poles. Consequently, with ten scattering channels one has 
for each eigenvalue of the confinement spectrum 1024 poles in 1024 Riemann sheets, out of which 
usually only one contains relevant poles in each ^/i interval in between the thresholds. Those 
can be identified with the known resonances, like the p pole in tttt scattering, or the K* pole in 
Ktt scattering. However, there may always be a pole in a nearby Riemann sheet just around 
the corner of one of the thresholds, which can be noticed in the partial-wave cross section. The 
study of such poles is an interesting subject by itself [18,44]. 

Once the four model parameters are adjusted to the experimental phase shifts and cross 
sections, the pole positions can be determined and compared with tables for meson spectroscopy. 
For the purpose of the present investigation, we shall focus next on the singularity structure of the 
S matrix for the lowest-lying poles in 5'-wave meson-meson scattering, employing the simplified 
model of Sec. 3 above. 

5 Scattering-matrix poles 

In the hypothetical case of very small couplings for the three-meson vertices, we obtain poles in 

the scattering matrix that lie close to the eigenvalues of the confinement spectrum. Let us denote 
by Ml and M2 the meson masses, and by the difference between the complex- energy pole 
of the scattering matrix and the energy eigenvalue E^f of the nearby state of the confinement 
spectrum. Using formula (5), we obtain 

AE ^ {R(s) - il{s)} \J^Nf ■ (8) 
We may distinguish two different cases: 

(1) En > M1 + M2 (above threshold), 

(2) En < Mi + M2 (below threshold). 

When the nearby state of the confinement spectrum is in the scattering continuum, then AE 
has a negative imaginary part and a real part, since both R{s) and I{s) of formula (8) are real, 
and I{s) is moreover positive. The resonance singularity of the scattering matrix corresponding 
to this situation is depicted in Fig. 3. 
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Notice that the resonance pole is in the lower half of the complex-energy plane (second 
Riemann sheet), as it should be. 



SJm(-v/i) 










En cut 




threshold 




= Ml + M2 


Jae 








• resonance position 



Figure 3: When the confinement state on the real ^/s axis is in the scattering continuum, then for small 
coupling (perturbative regime) the resonance pole moves into the lower half of the complex ^/s plane. 

When the nearby state of the confinement spectrum is below the scattering threshold, then 
Ai? has only a real part, since I{s) turns purely imaginary below threshold, whereas R{s) remains 
real. The bound-state singularity of the scattering matrix corresponding to this situation is 
depicted in Fig. 4. 

Note that the bound-state pole is on the real axis of the complex-energy plane, as it should 

be. 











LsRe(v^) 


bound state position 







En 
m 


threshold cut 








= Ml M2 



Figure 4: When the confinement state on the real ^/s axis is below the lowest scattering threshold, 
then the bound-state singularity comes out on the real ^/s axis. 

6 Threshold behavior 

Near the lowest threshold, as a function of the overall coupling constant, S'-wave poles behave 
very differently from P- and higher-wave poles. This can easily be understood from the effective- 
range expansion [45] at the pole position. There, the cotangent of the phase shift equals i. 
Hence, for S waves the next-to-lowest-order term in the expansion equals ik {k represents the 
linear momentum related to s and to the lowest threshold). For higher waves, on the other hand, 
the next-to- lowest-order term in the effective-range expansion is proportional to k'^. 
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Poles for P and higher waves behave in the complex k plane as indicated in Fig. 56. The two 
/c-plane poles meet at threshold (A; = 0). When the coupling constant of the model is increased, 
the poles move along the imaginary k axis. One pole moves towards negative imaginary A;, 
corresponding to a virtual bound state below threshold on the real ^/s axis, but in the wrong 
Riemann sheet. The other pole moves towards positive imaginary A;, corresponding to a real 
bound state. 



bound 
state 1 


Iin(A;) 

Re(A:) 






virtual ^ 
bound 
state 


resonance 
pole 

S wave 



bound 
state 1 


Im(A;) 


virtual ^ 
bound 
state 


resonance 
pole 

P wave 



(a) 



Figure 5: Variation of the positions of scattering-matrix poles as a function of hypothetical variations in 
the three-meson- vertex coupling, for S waves (a), and for P and higher waves ih). The arrows indicate 
increasing coupHng constant. 



For 5'-wave poles, the behavior is shown in Fig. 5a. The two A;-plane poles meet on the 
negative imaginary k axis. When the coupling constant of the model is slightly increased, both 
poles continue on the negative imaginary k axis, corresponding to two virtual bound states below 
threshold on the real ^/s axis. Upon further increasing the coupling constant of the model, one 
pole moves towards increasing negative imaginary k, thereby remaining a virtual bound state 
for all values of the coupling constant. The other pole moves towards positive imaginary k, 
eventually passing threshold [k = 0), thereby turning into a real bound state of the system of 
coupled meson-meson scattering channels. Hence, for a small range of hypothetical values of the 
coupling constant, there are two virtual bound states, one of which is very close to threshold. 
Such a pole certainly has noticeable influence on the scattering cross section. 

Although we are not aware of any experimental data that could conflrm the above-described 
threshold behavior of poles, we suspect this to be possible for atomic transitions in cavities. 
Unfortunately, it does not seem likely that in the near future similar processes can be studied 
for strong coupling. 
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7 The low-lying nonet of A^-wave poles 



The nonet of low-lying 5"- wave poles behave as described in Sec. (5), with respect to variations 
of the model's overall coupling constant. However, they do not stem from the confinement 
spectrum, but rather from the cavity. For small values of the coupling, such poles disappear 
into the continuum, i.e., they move towards negative imaginary infinity [20], and not towards an 
eigenstate of the confinement spectrum as in Fig. 3. 
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Figure 6: Hypothetical movement of the Kq(800) pole in Ktt (/ = 1/2) 5-wave scattering as a function 
of the coupling constant A (A is increased by steps of 0.01 unit, starting from 0.9). The two branches on 
the imaginary k axis (see Fig. 5) both result in poles on the real axis in the E = ^/s plane. However, in 
order to visualize their movement as a function of the coupling constant, wc give virtual bound states 
{Im{k) < 0) a small negative and bound states {Im{k) > 0) a small positive imaginary part. In the 
inset we show in more detail how the poles coming from the lower half of the complex E = y/s plane 
end up on the real axis when the coupling constant is increased from 1.2 to 1.25. 

In Fig. 6 we study the hypothetical movement of the ^0(800) pole in Xtt S-w&Ye scattering. 
The physical value of the couphng constant equals 0.75, which is not shown in Fig. 6. A figure 
for smaller values of the couphng constants can be found in Ref. [20]. The physical pole in 
Kn isodoublet S-wme scattering, related to experiment [17], comes out at 727 — i263 MeV in 
Ref. [19]. Here we concentrate on the threshold behavior of the hypothetical pole movements in 
the complex k and planes as a function of the coupling constant. Until they meet on the axis, 
which is for a value of the coupling constant slightly larger than 1.24, we have only depicted the 
right-hand branch. 

The pole corresponding to the one moving downwards along the imaginary k axis moves to the 
left on the real ^/s axis. The pole which moves upwards along the imaginary k axis initially moves 
towards threshold and then turns back, following the former pole, but in a different Riemann 
sheet. In the inset we clarify the motion of the (S-matrix singularities just before and just after 
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they represent virtual and real bound states. Notice that, since we took 0.14 GeV and 0.50 GeV 
for the pion and the kaon mass, respectively, we end up with a threshold at 0.64 GeV. 

It is interesting to note that in a recent work by Boglione and Pennington [46] a zero-width 
state is found below the Kn threshold in S-wme scattering, instead of the Kq{800) resonance. 
Here, we would obtain such a state for unphysical values of the coupling. 

In Fig. 7 we have depicted the movement of the ao(980) pole in S wave 1=1 KK scattering 
(threshold at 1.0 GeV) on the upwards-going branch. One observes a very similar behavior as in 
the case of Kt: scattering, but with two important differences, to be described next. 
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Figure 7: Pole movement of the ao(980) as a function of the coupling constant A for KK [J = 1) 5-wave 
scattering (A is increased by steps of 0.01 unit, starting from 0.02). As in the case of the i^Q(800) pole, 
the two branches on the imaginary k axis (see Fig. 5) both result in poles on the real axis in the E = ^/s 
plane. However, in this case we only study the upwards moving branch. Also here, in order to visualize 
their movement as a function of the coupling constant, we give virtual bound states {I'm{k) < 0) a small 
negative and bound states {Im{k) > 0) a small positive imaginary part. In the inset we show in more 
detail how the poles coming from the lower half of the complex E = y/s plane end up on the real axis 
when the coupling constant is increased from 0.50 to 0.51. 



The iro(800) poles meet on the real ^/s axis only 16 MeV below threshold (see Fig. 6), and 
for a value of the coupling constant which is well above the physical value of 0.75, whereas the 
ao(980) poles meet 238 MeV below threshold, when the coupling constant only equals slightly 
more than 0.5. At the physical value of the coupling constant, the ao(980) pole is a real bound 
state some 9 MeV below threshold. 

But there is yet another difference. Whereas the Kt: channel represents the lowest possible 
scattering threshold for the ^^0(800) system, KK is not the lowest channel for the ao(980). In 
a more complete description, at least all pseudoscalar meson-meson channels should be taken 
into account. One of these is the 7777 channel, which has a threshold well below KK. Conse- 
quently, upon including the 7777 channel in the model, the pole cannot remain on the real ^ 
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axis, but has to acquire an imaginary part in a similar way as shown in Fig. 3. In Ref. [19] we 
obtained a resonance-hke structure in the rjir cross section, representing the physical ao(980). 
The corresponding pole came out at 962 — i28 MeV. 

For the /o(980) system the situation is very similar to that of the ao(980). Assuming a pure 
ss quark content [47], we obtain for the variation of the corresponding pole in KK (/ = 0) 
S-wavc scattering a picture almost equal to the one shown in Fig. 7. However, only in lowest 
order the KK channel can be considered the lowest threshold for the /o(980) system. In reality, 
ss also couples to the nonstrange quark-antiquark isosinglet, namely via the KK channel, and 
hence to tttt [48]. Nevertheless, this higher-order coupling turns out to be rather weak, which 
implies that the resulting pole does not move far away from the KK bound state. In Ref. [19] 
we obtained a resonance-like structure in the tttt cross section representing the physical /o(980). 
The corresponding pole came out at 994 — i20 MeV. 

At lower energies, we found for the same cross section a pole which is the equivalent of the 
Kq{800) system, but now in tttt isoscalar S-wme scattering. This pole at 470 — i208 MeV can 
be associated with the a or /o(600) meson, since it has the same quantum numbers and lies in 
the ballpark of predicted pole positions in models for the a (a complete overview of a poles can 
be found in Ref. [49]). 

We do not find any other relevant poles in the energy region up to 1.0 GeV. 

8 Summary and conclusions 

It should be clear from the foregoing that the light scalar mesons do allow a description on the 
basis of normal qq states, provided one accounts for mesonic loops. The crucial point is that, 
while for the other mesons the effect of coupling to open and closed meson-meson channels is 
more or less perturbative, giving rise to complex or real mass shifts of the qq confinement states, 
in the case of the scalar mesons the very strong ^Pq coupling to S-wme pseudoscalar-pseudoscalar 
channels forces additional, highly non-perturbative poles in the S matrix to approach the physical 
region, unhampered by any centrifugal barrier in the two-meson sector. In other words, these 
singularities do not stem from the confinement spectrum, at least not in a simple way, but 
are rather due to the ^Pq barrier providing the communication between the qq states and the 
meson-meson continuum. Nonetheless, these poles give rise to very pronounced structures in the 
5"- wave scattering of pseudoscalar mesons, which can be either clear-cut resonances (ao(980) and 
/o(980)), or broad non-resonant bumps (/o(600) and Kq{800)). 

These conclusions we were able to make more quantitative by employing a simplified two- 
channel model with a delta-shell transition potential, allowing the derivation of a general, closed- 



15 



form expression of the scattering phase shift for an arbitrary confinement mechanism. This for- 
mula, which we baptized "Resonance-Spectrum Expansion" , turns into a standard Breit-Wigner 
form in the vicinity of a particular resonance, in the limit of small coupling and width, whereas it 
describes hard-sphere scattering in the limit of large coupling. Application of this formalism to 
the specific case of the scalar mesons unmistakably shows that all states of the light nonet have 
the same origin, and, in principle, the same threshold behavior. It just happens that the various 
thresholds of the different scalars are very disparate, not in the least due to the small pion mass 
and hence to chiral symmetry. 

In particular, we found, in the complex-energy plane, a nonet of S'-matrix poles representing 
the ao(980), /o(980), /o(600), and K*{800), the latter two having real parts of 0.47 GeV and 
0.73 GeV, respectively, and imaginary parts of 0.21 GeV resp. 0.26 GeV. Whether or not these 
poles manifest themselves as clear physical resonances [50] is not so relevant here in view of their 
common nature, as we have demonstrated in detail for the Kq{800) and ao(980). Besides these 
"non-perturbative" states, we also found the confinement-ground-state nonet of scalar mesons 
/o(1370), ao(1450), -^"0(1430), and /o(1500). The latter poles vary as a function of the coupling 
constant exactly the way indicated in Fig. 3. For vanishing coupling, they end up on the real 
energy axis at the positions of the various ground-state eigenvalues of the confinement spectrum, 
which are the light-flavor ^Pq states at 1.3 to 1.5 GeV [51-53]. Recall that, in contrast, the poles 
of the light scalars move towards negative imaginary infinity in this limit. 

In conclusion, we should mention that a coupled-channel or unitarized approach to the scalar 
mesons, similar in spirit to our present and previous [19, 34] works, is rapidly gaining advocates 
[54-57]. Nevertheless, in none of these works a simultaneous discription of the two lowest scalar 
nonets is achieved as outhned above, most notably the still widely contested [50,54,55] but now 
experimentally confirmed i^^^(800). As a final remark, let us point out that in our full model with 
many coupled channels [19], all channels contribute to the states under the resonance, and not 
just one specific channel. However, both the full and the simplified model produce very similar 
masses and widths for the light scalars when reproducing the experimental phase shifts in the 
relevant energy region. This lends additional quantitative support to our predictions for these 
observables presented here. 
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A The T matrix for meson-meson scattering 



In Eq. (1), we must eliminate ^c, since it vanishes at large distances and is thus unobservable. 
Formally, this can be done in a straightforward way. We then obtain the relation 



(9) 



By comparison of Eq. (9) with the usual expressions for the scattering wave equations, we 
conclude that the generalized potential V is here given by 

V = [Vtf {E-H,)-' Vt . (10) 
The matrix elements of the T-operator are defined by the Lippmann-Schwinger equation 



t(p,p' •,z) = V(p,p') + 
+ J d^k' J d^kv[p,k^ Gf (k\k]z^ V{k ,p') + 

+ J d^k'" J d^k" J d^k' j d^k V [p,k" ) Gf (%"' , k" ; V {%" , k^ X 

xGf(k\k-z^ V{k,p') + ... , (11) 

where the Green's operator Gf{z) corresponds to the self-adjoint free Hamiltonian Hf, according 
to 



Gf Ik ,k;z 



k 



>-Hf)-' 



k 



'-I'J 



k 



k 



2^fZ — k'^ 

Substitution of relation (12) in expression (11) yields for the T-matrix elements the form 
t(p,p';z) = V{p,p') + Jd'kV{p,k) ^-^^^V{k,p') + 



(12) 



+ 



d^kV{p,k 



2f,f 



V k ,k 



(13) 



V{k,p') + 



2fifZ-k'^ \ ' / 2i2fZ-P 
At this stage, we must make a choice for the operators He, describing the confinement dy- 
namics in the interaction region, Hf, representing the dynamics of the scattered particles at large 
distances, and Vt, which stands for the transitions between these two sectors. For an arbitrary 
spherically symmetric confinement potential Vc, we define these operators in configuration space 

by 
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He = 



2iic 



+ niq + rriq + Vc{r) , Hf = 



V2 

+ Ml + M2 



and Vf 



A 



2iica 



S (r — a) 



(14) 



The various mass parameters of Eq. (14) are defined in Table 3. Tfie transition potential Vf, 
which provides the communication between the confined channel and the scattering channel, is 
an extreme simplification of potentials [42] that may describe the breaking of the color string. 
Here it is assumed to only act when the particles are at a distance r — a, thus having the form 
of a sperical delta shell. 



symbol 


definition 


/^/ 


constituent (anti-) quark mass 
reduced mass in confinement channel 
meson masses 

reduced mass in scattering channel 



Table 3: Definition of mass parameters used in Eq. (14). 

In configuration space, we may then write the non-relativistic 2x2 stationary matrix wave 
equation (1) in the form 



A 



2/ica 

A 
2//ca 



5{r-a) ipf{f) , 

6 {r — a) ipc{^) 



A.l The Born term 

In the momentum representation, Eq. (10) takes the form 



V{p,p') = {p\v\p') = {p\ [Vtf {E{p)-H,r' Vt 



P 



(15) 



(16) 



The total center-of-mass energy E and the linear momentum p are, through Eq. (14), related by 



Eip) 



P 



+ Ml + M2 



(17) 
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We denote the properly normalized eigensolutions of the operator He (14), corresponding to the 
energy eigenvalue E„i, by 



(f \n£m) = ll^) (f) J^Ur) , with n ^ 0, 1, 2, ^ 0, 1, 2, m ^ -£,...,+£ . (18) 
So, by letting the self-adjoint operator He act to the left in Eq. (16), we write 



{p\V 



P 



^ (p I [Vtf \nim) {nem\ {E{p) - H,)-' Vt 

nim 



\nim) {nim\ 



Next, we insert several times unity to obtain 



P 



(19) 



(p^l V \p') = Jd^^ J d^'r" J Sr'" j Sr' 

nlm 



1 



(p|f)(f| [Vt] f'j{j''\nlm){nlm 



-•III \ / -till 



f')(f' 



P 



E{p) - Er^e 

Two of the four integrations are trivial, since the local transition potential has the form 

A 



(f I Vt r') = -— 5{r-a) S^^^ (f - f ') 



(20) 



(21) 



By inserting expression (21) into Eq. (20), also substituting {f \p) = e^P ' ^/(27r)^/^, we get 



{p\ V \p') = j j d^^' 
1 e-'^P-^ A 



(22) 



X 



5 {r — a) {f \nim) (nim 



■ ->i -ti 

^/\ A e^P -r 

r' ) 5 ir' — a) 



2/Xca 



(27r)=^/2 



E{p) - Er,e (27r)3/2 2fi,a 
Next, we observe that the radial parts of the two remaining integrations are also trivial, because 
of the two delta functions. So we twice insert the expression for the confinement eigenfunctions 
of Eq. (18), to obtain 



{p\V 



P 



A 



(27r)^ \2iica^ 



J2 JdQa^ J dfl' 

nim 



(23) 



X 



E{p) - E„e 



e-^^-«^rW(f)^„,(a)rF(r')^:.(a)e^^ • 



ar 
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For the integrations over the angles we introduce Bauer's formula [58] , 



-ik ■ f 



resulting in 



/ 



^-ip ■ ar ^ 



> = 471 je{pa) 



HYy^Hv) 



Substitution of the relations (25) into Eq. (23) leads to the expression 



p 



AV ^ 

2^ }tra 



1 



(24) 



(25) 



Eip) - E^e 



Y^^* (p) Y^^ (p') U{pa) jeip'a) |^„,(a)|^ , (26) 



where the summation over the magnetic quantum number m can be performed by the use of the 
addition theorem, thus shaping the Born term (16) into its final form 



{p\V 



2^2 oo 



J2{2i + 1) Pe {p-p') J, {pa) jeip'a) 



"o E{p) - E^i 



(27) 



A. 2 The second- order term 



For the second-order term, we start by substituting the result (27) into the second term of 
expansion (13), giving rise to the expression 



J ol^ fJ'c i=Q n=0 ^\r) ^n. 



|2 



The details of the k integration are discussed in Appendix (A. 5). Since, as required by Eq. (34) 
below, E{k) is quadratic in k, we find for expression (28) the result 
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X 



4 



|2 



oo 



A. 3 To all orders 

Following steps similar to those in Appendix (A. 2), it is now straightforward to determine the 
higher-order contributions to the expansion (13). For the full T matrix to all orders, one ends 
up with the result 

t(p,p') = (30) 



A. 4 Scattering matrix and phase shift 

For radially symmetric interactions, it is useful to define the partial-wave matrix element of 
the on-shell (p ' = p) T matrix, according to the relation 

oo 

T{p) = E(2£+1)P,(P-P') Te{p) . (31) 

Hence, also using the result of Eq. (30), we find for the partial- wave scattering amplitude Si{p) 
the expression 



Seip) 



8i7r^l^fpTi{p) 



(32) 



>M^fP 



Xa 



ji{pa) hf\pa) J2 



^0 ^iP) - ^ne 
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For the partial- wave scattering phase shift defined by S^ip) — , one obtains 

from Eq. (32) the result 



cotg(5.(p)) = — r° ,1 .2 ■ (33) 



Formula (3) for the meson-meson scattering phase shift is based on the latter equation, but 
not exclusively, as we discuss in Appendix (C). 

A. 5 Details of the k integration 

Let us study the momentum-space integration 

{p,P ;^^-Je) = (34) 

= 2^ / dn, p,(p- k ) p, (k.p') f edk ''^I'^J^fh , [k-) . 

where fi represents an arbitrary well-behaved function of fc^. For the integration over the angles, 
we can employ the orthogonality of spherical harmonics. Hence, we must concentrate on the 
radial integration, i.e.. 

We shall show below that the integration can easily be performed, yielding 

by using the following properties of the spherical Bessel and Hankel functions: 

ji(e^^ ka) ^ e''''^ ji{ka) and /ij'^ (e^^ A;a) = e"^^^ /if (A;a) . (37) 

For large imaginary part of the argument ka, the function h!l\ka) tends to zero. Therefore, we 
can close the integration path in the complex k plane by a non-contributing semicircle in the 
upper half plane. If we then set 2^z — {p + ief, taking the limit e J, after the integration, the 
integral (36) can be simply computed with Cauchy's residue theorem, yielding 
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2 n{po) h^Pipa) fi(f) . 
Putting everything together, we obtain for Eq. (34) the final result 



(38) 



(39) 



B The phase shift in the configuration-space approach 



The radial wave equation, following from Eq. (15) by choosing 



is given by 



ijjf (r ) 

he XV(r) 
XV{r) hf 



(40) 



Uc 



= E 



Uc 



(41) 



where (see Eq. (14)) 



he 



hi 



V 



2 Ate V '^^^ 
1 / (P 



'ifif \ dr 
1 



+ 



+ 



+ 1)' 



+ ruq + ruq + Vc{r) , 



+ Ml + M2 , and 



2/Xca 



5 {r — a) 



(42) 



In this section, we study the solutions of the 2x2 radial wave equation (41) in configuration 
space. 

For r < a as well as for r > a, one has to solve the two uncoupled differential equations, 
equivalent to A = 0, given by 



he \ / Ue 
hf ) [uf 
At r — a one has the boundary conditions 



E 



Ue 



(43) 
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and 



1 /_ d Ucjr) 
\ dr 



+ 



d uJr) 



ria 



dr 



A , 1 (_ d Uf{r) 



(44) 



+ 



(i Uf{7 

dr 



(45) 



Uc{r T a) = iic('^ i a) 

Uf{r t 0.) = %('^ i a) 

Further boundary conditions are the usual ones: the wave functions and uj must both vanish 
at the origin. Towards infinity, the wave function Uc must be damped exponentially, since Uc 
describes a confined system, whereas, above threshold {E > Mi + M2), the wave function Uf 
must have an oscillating behavior, describing the scattered mesons. 

Let us denote by Fc^^ and Gc,e the solutions of the upper differential equation of Eq. (43) , which 
vanish at the origin and fall off exponentially at infinity, respectively. The solution Fc^^{E, r) is 
regular at the origin, but is irregular at infinity, except for some special cases known as the 
confinement spectrum, whereas Gc,e{E, r) behaves regularly at infinity but not at the origin. For 
arguments belonging to the confinement spectrum, F^/ and Gc/ are degenerate and normalizable. 

The lower differential equation Eq. (43) is solved by the functions J and N defined by 



Ji{k,r) = k ^rj^{kr) and N({k,r) = k^^^rni{kr) , 
where the linear momentum k is given by 

k'' = 2iif {E - Ml - M2) , 



(46) 



(47) 



and where ji and ng are the spherical Bessel and Neumann functions, respectively. The solution 
Jii^k, r) satisfies the usual boundary condition at the origin, i.e., J^(r — > 0) ^ 0. For the solutions 
(46) one has the Wronskian relation 



W{Mk,a),Ni{k,a)) = 



dNi{k,r) dJe{k,r) 
Je{k, r) Ni{k, r) 



= 1 



(48) 



dr dr 

A general solution of the Schrodinger equation (42), which, moreover, satisfies the boundary 
conditions at the origin and at infinity, reads 
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UciE,r) 
Uf{E,r) 



Je{k,r) Af 



r < a 



(49) 



r > a 



Je {k, r) A;2^+icotg {Si{E)) - TV, (A;, r)] Bf ^ 

where Ac, Af, Be and Bf are normalization constants, which are not independent because of the 
boundary conditions (44) and (45). From the boundary conditions (45) we derive the equations 



F,,e{E,a)A, = G,,e{E,a) B, , 



J, {k, a) A;2^+icotg {6e{E)) - N, {k, a) 



Je {k,a) Af = 
Similarly, from the boundary conditions (44) we obtain 

^JG[iE,a)B,-F^iE,a)A,: 



Bf 



(50) 



2^Ji{k,a)Af , 



^([j',{k,a)k''^'cotg{6e{E))-N',{k,a)]Bf-J^{k,a)Af) = ^F,,e{E,a)A, 



(51) 



By elimination of the normalization constants, one finds for the cotangent of the phase shift the 
expression 



k''+'cotg{Se{E)) = 



Ne {k, a) 



Je {k, a) 

Using moreover the relations (46), we arrive at 



n -1 



W{F,^,,G,,e) 



(52) 



cotg(5,(E)) 



{ka) 



f,. F,,,{E,a)G,,,{E,a) 



^2 l^fk ^2 



-1 



(53) 



je (ka) 

Note that, in Eq. (53), we have obtained a simpler expression for the phase shift than in 
Eq. (33), in particular numerically, since the convergence of the sum in Eq. (33) is very slow. It 
is, moreover, amusing that we have thus derived, as a side result, a "simple" (depending on one's 
taste) proof for the following spectral representation of the Green's function: 



E 

n=0 



2//e F,,e{E, a)G,,e{E, a) 



E-E„e WiF,^,iE,a),G,,eiE,a)) ' 



(54) 



This relation must be true for any confining potential with an infinite discrete set of radial and 
angular excitations. The set of functions {J-'n£ n = 0, 1,2, . . .} represents, for orbital angular 
momentum £, a full set of radial eigensolutions, with eigenvalue Ene, of the Hamiltonian He- 
Furthermore, Fc^£ and Gc,e represent two linearly independent solutions for any value of the 
energy E. 
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C More realistic "phase" transitions 

The results of Appendices (A) and (B) are illustrative for studying the skeleton of the calculus 
involved in determining the elastic scattering phase shifts, but do not represent the coupling of 
quark-antiquark states to meson-meson systems. An easy way to understand this is by comparing 
the parities of both phases of the system. Due to Fermi statistics, the parity of the qq phase is 
given by P = (—1)^=+^. On the other hand, from Bose statistics follows that the parity of the 
meson-meson phase is given by P = (—1)^^. One concludes that the orbital angular momenta Ic 
and for respectively the qq phase and the meson-meson phase, must differ at least one unit. 
Hence, the potential given in Eq. (14) cannot couple the two different phases. Nevertheless, 
expression (33) for the scattering phase shift is a powerful tool for further investigation. In 
Appendix (B) we obtained a different form for this quantity. Now the procedure leading from 
Eq. (41) to Eq. (53) can be repeated for the more general case, given by 



K = 7^ ( -33 + — — ^2 + rriq + ruq + V;(r) , 



1 / (f £/(£/ + l)\ 

V = ^M-a) . (55) 

Substituting £ in Eq. (53) conveniently by either £c or if , one gets for the scattering phase shift 
the result 



+ fx ^^^^ (ka) 



jif (ka) 

Then we can make use of Eq. (54) to arrive at 



(56) 



cotg [Si^{E)) = . . - A -^-^j {ka) ^- — - 

^ ' ^" '''" ^ I ^l^c n=0 E ~ ^nl^ 



(57) 



-.2 00 I -77 ( „\ |2" 

jif (ka) 

This is the result that can be used for our model. However, note that we have not specified how 
we arrived from Eq. (15) to the relations (55). Operators which provide for the communication 
between channels of different orbital angular momenta and quark-antiquark spins, can be con- 
structed [23,42,59], but the procedure involves rather technical recoupling schemes of quantum 
numbers. We do not intend to go into the details here. 

Extension of the formahsm to potentials more realistic than a spherical delta shell is also 
straightforward, and can be found in Refs. [35,60]. However, it is opportune to mention here 
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that the solutions of the full model [19,23] do not behave differently from the simple formula (57) 
in the energy domain under study in this work, only yielding some minor differences in the 
numerical results. 
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